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Abstract 

We show that radiation pulses from an X-ray Free-Electron Laser (XFEL) with a 
planar undulator, which are mainly polarized in the horizontal direction, exhibit a 
suppression of the vertical polarization component of the power at least by a factor 
A^/(47rLg)^, where A^, is the length of the undulator period and Eg is the FEL field 
gain length. We illustrate this fact by examining the XFEL operation under the 
steady state assumption. In our calculations we considered only resonance terms: 
in fact, non resonance terms are suppressed by a factor A^/(47rLg)^ and can be 
neglected. While finding a situation for making quantitative comparison between 
analytical and experimental results may not be straightforward, the qualitative 
aspects of the suppression of the vertical polarization rate at XEELs should be easy 
to observe. We remark that our exact results can potentially be useful to developers 
of new generation EEL codes for cross-checking their results. 
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1 Introduction 


A Free-Electron Laser (FEL) amplifier which starts up from the shot noise in 
the electron beam is known as a self-amplified spontaneous emission (SASE) 
FEL. In the SASE case, the amplification process has its origins in the den¬ 
sity finctnations of the electron beam. SASE FELs are capable of prodncing 
coherent, tnnable FEL radiation down to a fraction of an Angstrom BIS. A 
SASE FEL which operates in the X-ray wavelength range is called XEEL. 

With the advent of EELs, X-ray radiation pnlses with nnprecedented charac¬ 
teristics were made available to the scientific commnnity. Compared to conven¬ 
tional synchrotron radiation sonrces. X-ray EELs (XEELs) offer an increase 
in peak brightness of many orders of magnitnde as well as nltrashort pnlses 
in the femtosecond time scale. In this paper we consider an additional featnre 
of XFEL pnlses, which is nsefnl in many experiments. Namely the fact that 
the pnlses prodnced by an XFEL with horizontal planar nndnlator exhibit an 
extremely small component of the electric field in the vertical direction. In 
particnlar, here we will show that for a typical XEEL setnp the horizontally 
polarized component of radiation is greatly dominant, and that only less that 
one part in a million of the total intensity is polarized in the vertical plane. 

The stndy of XFEL polarization characteristics is obviously deeply related 
to the problem of electromagnetic wave amplification in XEEL, which refers 
to a particnlar class of self-consistent problems. It can be separated into two 
parts: the solntion of the dynamical problem, i.e. finding the motion of the 
electrons in the beam nnder the action of given electromagnetic fields, and 
the solntion of the electrodynamic problem, i.e. finding the electromagnetic 
fields generated by a given contribntion of charge and currents. The problem 
is closed by simnltaneons solntion of the field eqnations and of the eqnations 
of motion. 

Let ns consider the electrodynamic problem more in detail. The eqnation for 
the electric field follows from Maxwell’s eqnations. One obtains, in Gaussian 
nnits: 



( 1 ) 


With the help of the identity 


V X (V X E) = V(V ■ E) - 


( 2 ) 
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and Poisson equation 


V ■ -E = dvrp 


( 3 ) 


we obtain the inhomogeneous wave equation for E 



4:71 p + 47r— . 


( 4 ) 


Once the charge and current densities p and j are specified as a function of 
time and position, this equation allows one to calculate the electric field E at 
each point of space and time |2|. The current density source provides the main 
contribution to the radiation held in an PEL ampliher, and the contribution 
of the charge density source to the amplihcation process is negligibly small. 
This fact is commonly known and accepted in the FEL community. However, 
we have been unable to hnd a proof of this fact in literature, except book |3] 
and review [3], which are only the publications, to the authors’ knowledge, 
dealing with this issue. 

Due to linearity, without the gradient term the solution of Eq. (4) exhibits 
the property that the radiation held E points in the same direction of the 
current density j. An important limitation of such approximation arises when 
we need to quantify the linear vertical held generated in the case of an XFEL 
with planar undulator. In the case j points in the horizontal direction (for a 
horizontal planar undulator), according to Eq. (4), which is exact, only the 
charge term is responsible for a vertically polarized component of the held: if 
it is neglected, one cannot quantify the linear vertical held anymore. 

Similar to the process of harmonic generation, the process of generation of the 
vertically polarized held component can be considered as a purely electrody¬ 
namic one. In fact, the vertically polarized held component is driven by the 
charge source, but the bunching contribution due to the interaction of the elec¬ 
tron beam with the radiation generated by such source can be neglected. This 
leads to important simplihcations. In fact, in order to perform calculations of 
the radiation including the vertically polarization component one can proceed 
hrst by solving the self-consistent problem with the current source only. This 
can either be done in an approximated way using an analytical model for the 
FEL process or, more thoroughly, exploiting any existing FEL code. Subse¬ 
quently, the solution to the self-consistent problem can be used to calculate 
the hrst harmonic contents of the electron beam density distribution. These 
contents enter as known sources in our electrodynamic process, that is Eq. 
(4). Solution of that equation accounting for these sources gives the desired 
polarization characteristics. 
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Analytical descriptions allow for a proper understanding of the physical prin¬ 
ciples underlying the phenomena under study, and also provide convenient 
testing of numerical simulation codes. A SASE XFEL is rather difficult to 
be described in a fully analytical way. In all generality, its radiation can be 
represented as a non-stationary random process, and its analytical description 
is further complicated by the fact that the electron bunch combines both the 
features of the input signal and of the active medium with time-dependent 
parameters. 


Approximations particularly advantageous for our theoretical analysis include 
the modeling of the electron beam density as uniform, and the introduction of 
a monochromatic seed signal. Realistic conditions satisfying these assumptions 
are the use of a sufficiently long electron bunch with a longitudinal stepped 
prohle and the application of a scheme in the SASE mode of operation for 
narrowing down the radiation bandwidth. In the framework of this model it 
becomes possible to describe analytically all the polarization properties of the 
radiation from an XFEL H 


The simplicity of our model offers the opportunity for an almost completely 
analytical description of an XFEL in the linear mode. As remarked above, 
a complete description of the operation of an XFEL can be performed only 
with time-dependent numerical simulation codes. Application of the numer¬ 
ical calculations allows one to describe the most general situation, including 
arbitrary electron beam quality and nonlinear effects. Finding an analytical 
solution is always fruitful for testing numerical simulation codes. Up to now, 
in conventional EEL codes, the contribution of the the charge source assumed 
to be negligible small. However, the charge term alone is responsible for the 
vertically polarized radiation component, which is our subject of interest. Our 
analytical results for the high-gain linear regime are expected to serve as a 
primary standard for testing future EEL codes upgrades. 


^ Our model can be close to real situations, based on two techniques recently real¬ 
ized at the LCLS. The first is a technique for producing a uniform electron bunch, 
and is heavily relying on the use of a slotted spoiler foil in the last bunch compressor 
chicane [6]. The method takes advantages of the high sensitivity of the FEL gain 
process to the transverse emittance of the electron bunch. By spoiling the emit- 
tance of most of the long nonuniform bunch while leaving short unspoiled temporal 
slice, one can produce electron bunch with relatively uniform active medium. The 
second technique is self-seeding, an active filtering technique allowing to narrow the 
FEL bandwidth down to almost the Fourier limit [7j. Combination of these two 
techniques is now widely used at the LCLS. 
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2 XFEL radiation in resonance approximation 


As has been discussed in the Introduction, a quantihcation of the linear ver¬ 
tical held expected in FELs can be made on the basis of an electrodynamical 
model. We will restrict our attention to FELs driven by planar undulators. A 
more extended treatment of higher harmonic numbers is conceptually identical 
to the treatment in this paper, and only differs as regards actual calculations. 
Proper initial conditions following e.g. from start-to-end simulations are given 
as input to an FEL self-consistent code, which calculates the electron beam 
bunching from the interaction of the beam with the hrst harmonic radiation. 
As discussed above, in FEL codes the hrst harmonic is conventionally calcu¬ 
lated accounting only for the leading terms under the resonant approxima¬ 
tion, which yields perfectly linearly polarized radiation. Therefore, the output 
of these codes cannot be directly used to quantify the linear vertical held of 
FELs. However, the electron beam bunching can be inserted into the held 
equation as sources. A relatively simple electrodynamical model can then be 
developed in order to calculate corrections to the leading resonant terms, and 
therefore to calculate the linear vertical held. 

Paraxial Maxwell’s equations in the space-frequency domain can be used to 
describe radiation from ultra-relativistic electrons (see e.g. m- Let us dehne 
temporal Fourier transform pairs as: 


OO 

f{u) = J dt f{t) exp{iut) ^ f{t) 

— OO 



duf{oj) exp{—iojt) 


(5) 


and let us call the transverse electric held in the space-frequency domain, i.e. 
the Fourier transform of the real electric held in the time domain, E±{z, 
where f± = xe^ + yey identihes a point on a transverse plane at longitudinal 
position ; 2 , Cx and Cy being unit vectors in the transverse x and y directions. 
Here the frequency u is related to the wavelength A by cj = ‘I'xcjX, c being 
the speed of light in vacuum. From the paraxial approximation follows that 

the electric held envelope E± = E± exp [—iuz/c] does not vary much along 2 ; 
on the scale of the reduced wavelength A/(27r). As a result, the following held 
equation holds: 


V 


E_i{z,f±,u) 


g(z,f±,u) , 


( 6 ) 


where the diherential operator V is dehned by 
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(7) 


V = 


v,2 


+ 


2iuj d 


dz 


V_L^ being the Laplacian operator over transverse cartesian coordinates. Eq. 
(6) is Maxwell’s equation in paraxial approximation. The source-term vector 
g{.z, r_L) is specihed by the trajectory of the source electrons, and can be written 
in terms of the Fourier transform of the transverse current density, rj_, ca), 
and of the charge density, p[z,r±,uj), as 




g = -47r —j^ - V±p exp 


iuiz 

c 


( 8 ) 


In this paper we will treat jj_ and p as macroscopic quantities, without in¬ 
vestigating individual electron contributions, j^ and p are regarded as given 
data, that can be obtained from any FEL code. Codes actually provide the 
charge density of the modulated electron beam in the time domain p{z^ fr, ^)- 
A post-processor can then be used in order to perform the Fourier transform 
of p, that can always be presented as 


p = p{z, 


r;_L(z),a;)exp 


iuj 


So{,z) 

Vo 


(9) 


where ro±{z), So{z) and Vo are the transverse position, the curvilinear abscissa 
and the velocity of a reference electron with nominal Lorentz factor jo that is 
injected on axis with no deflection and is guided by the planar undulator held. 
Such electron follows a trajectory specihed by fo±{z) = Tqx^x + VoyVy with 


rox{,z) 


K 


cos{ku,z) 


r„y(z) = 0 . 


( 10 ) 


The corresponding velocity is described by Vo± {z) = Vox^x + VoyCy with 


Vox{,z) = -sin(/c^z) , Voy{z) = 0. (11) 

7o 

In Eq. (10), K = X^eHyj/{2TimeC^) is the undulator parameter, = 27r/fc^ 
being the undulator period, (—e) the negative electron charg^^ H^, the max¬ 
imal modulus of the undulator magnetic held on-axis, and rrie the rest mass of 
the electron. In writing Eq. (9) we are assuming that p varies approximately 
as exp [iujSo{z)/v^ over a distance z of order of several undulator periods. 

^ The minus sign on the right hand side of Eq. (9) is introduced for notational 
convenience and in agreement with the minus sign of the electron charge. 
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We note that for a generic motion one has 


u 


'So{z 2 ) - So{zi) Z2-Z1 


Z2 


= dz 


CO 


zi 


2072 ( 2 ;) 


( 12 ) 


where jziz) = 1/— Voz{zy/c^ and Voz{z) = \Jv‘l — Vo±{zY. Thus, with the 
help of Eq. (9), Eq. ( 8 ) can be presented as 


g = —dTrexp 


2 



dz 


u 

272 ( 7)0 


iu ^ 

-^Vo± 

& 



p{z,f^ - Toliz)^^) 


1 


(13) 


where we used the fact that = Vo^p- In fact, for each particle in the beam the 
relative deviation of the particles energy from is small, i.e. d'yj'yo ^ 1 - 

Therefore we can neglect differences between the average transverse velocity 
of electrons (hj_) and Eox, so that j_|_ = (hl)p ~ Vo 2 _p. Because of this we will 
also drop the subscript ” 0 ” in 70 . 


With the help of Eq. (13), Eq. ( 6 ) can also be presented as: 




—dTT exp 
where we have set 


i J dz 

L 0 


u 


-l 7 Vo± - V_L 


dz 


u 


2072(2;) I Lo 2 


uz uK^ 


p{z,f±-ro±{z),uj) , (14) 


2cj^(z) 472 — 0 


■ sin {2kyjz). 


(15) 


With the aid of the appropriate Green’s function an exact solution of Eq. 
(14) can be found without any extra assumption about the parameters of the 
problem: 


E±{z^f±^ijj) = / dz'- 


z — z' 


dr' 


lOJ ^ , ,, 


xp( 2 ;', r'j_ — ro_L( 2 :'), co) exp < ioj 


1 12 

r_L — r_L I 

2 c(x - z') 


+ i 


dz 


u 


2072 (^) 



( 16 ) 


where V'j_ represents the gradient operator with respect to the source point, 
while indicates the observation point. Integration by parts of the gra¬ 

dient terms leads to 


E_^{z,f±,uj) = — / dz '—^ - , 

c J z — z’ J V 

—oo ^ 

xp{z',r'i_ - roL{z'),uj) exp \i^T{z\r'^,uj) 


) ( Vo±{z') _ f± - r'j 

c z — z' 


(17) 


where the total phase 'hr is given by 


•hT — 


dz 


u 


2c-fl{z) 


+ U) 


\ t ± — r'_Lp 
2c{z - z') 


(18) 


We now make use of a new integration variable I = r'± — fo±{z') so that 


~ ^ iu f , 1 

E±{z,r±,u) = — / dz - - 

c J z — z' 

— oo 

xp(2;',/, w) exp \i^T{z\Puj 


dl 


Vo±{z') r^-fo±{z')-l 


z — z 


(19) 


and 


$7^ — 


dz 


UJ 


2c'^l{z) 


+ OJ 


\r± - rp^jz') - /p 
2c{z - z') 


( 20 ) 


Introducing the far zone approximation and making the appropriate substitu¬ 
tions into Eq. (19) yields the following held contribution calculated along the 
undulator; 


lU 


Ej_( 2 ;, rj_, w) =- dl dz'p{z',l,u)exp i^xiz',l,uj) 


cz 


X 


'K 


— sin {kujz') + 9x + QyCy 


7 


where 


( 21 ) 


9 



$7^ — ^ 






sin {2kwz') 


K9, 

iKc 


cos {k^z') 


+UJ I 7-^- {Oxlx + Oyly) + [6l + 9l) — 


kwjc c 


2c 


( 22 ) 


Here 9^ and 9y indicate the observation angles xjz and yjz. Moreover, since 
in Eq. (21) we introduced explicitly the trajectory inside the undulator, we 
need to limit the integration in dz' to a proper range within the undulator. 
We assume that this is done by introducing a proper function of z' as a factor 
p, which becomes zero outside properly dehned range, thus effectively limiting 
the integration range in z'. 


In this article we are interested in considering helds and electromagnetic 
sources originating from an FEL process. Imposing resonance condition be¬ 
tween electric held and reference particle, the self-consistent FEL process au¬ 
tomatically restricts the amplihcation of radiation at frequencies around the 
hrst harmonic: 


0J\q 2/c^cy^ , 


(23) 


where 


7 


2 

2 : 


1 + iLV2 ’ 


(24) 


and at emission angles 9 l^^^ l/b^- Our focus onto FEL emission also ex¬ 

plains the dehnition in Eq. (9). In fact, introduction of p is useful when p is 
a slowly varying function of on the wavelength scale. If the charge density 
distribution under study originates from an FEL process a stronger condition 
is satished, namely p is slowly varying on the scale of the undulator period 
and, as the FEL pulse itself, is peaked around the fundamental Uio- The 
words ’’peaked” or ’’around” the fundamental mean that the bandwidth is 
Aca/cuio 1. We quantify ”how near” the frequency u is to uio introducing 
the detuning parameter C, dehned as 


C 


u 

2c7f 


k'li 


Aoo , 

- k, 


(25) 


with Auj = uj — ujio. The detuning parameter C should indeed be considered as 
a function oi z, C = C(z) because, in general, the parameters of the undulator 
my be a function of ; 2 . All other dependencies on z, for example due to the fact 
that the energy of particles actually deviates from 7 and actually decelerate 
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during the FEL process, is accounted for in p. Indicating with we seek 
to calculate the hrst harmonic contribution at frequencies u around uio by 
making use of the well-known expansion 

CXD 

exp [iasin ('?/’)] = E Jp{a) exp [ip^j] , (26) 

p=—oo 

where Jp indicates the Bessel function of the hrst kind of order n. This allows 
to cast Eq. (21) into a form that is still fully general, but more convenient 
for enforcing the conditions, discussed above, related with FEL emission. One 
obtains 


OO OO OO 


Ei = 


cz 


dlx / dly / dz'p{z'{z\uj) eyip[i^c 


— OO —OO —OO 

OO OO 

X E E Jm{u)Jn{v)fiyZp 


m=—OO n=—OO 


imi 


X ■ 


K f 

( exp {i[R^ + l]k^z'} - exp {i[R^ - l]k^z'} 


+6x exp {iR^ky^z'} 




6y exp {iRyjkyjZ'} 


-'V f ’ 


(27) 


where 


ca 


Ruj = - n — 2m , 

cai 


(28) 


with 


1 


COlc 


1 + 7z (^x + 


(29) 


Moreover 


u = 


[1 - KV(4j^)] 




ui 4 


1 + 7^72 + 72(02 + 02)] 2(2 + 772 )’ 


V = 


c+o 2777 [1-777(472)] 0, 47770, 


1 + 772/2 + 72 (02 + 02 


2 + 772 


(30) 

(31) 


and 


4*0 = ^Ic 




(32) 
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Since the FEL process imposes that the magnitude of the charge density p 
does not vary much in z' over a period of the undulator A^, it follows that the 
fast oscillations in the exponential function in the integrand of Eq. (27) tend 
to suppress the integral unless = 0 , = — 1 or = 1 , that is when at 

least one of the exponential function is simply unity. Using Eq. (25) and Eq. 
(29)we rewrite Eq. (28) as 


Ru, = 


-r 

^11) 


l + 7f 




+ e: 


— n — 2m 


(33) 


Invoking again the FEL process, inspection of Eq. (33) shows that in the limit 
for C kw and <^1/71 the conditions above correspond to n = 1 — 2 m, 
n = 2—2m and n = —2m respectively. Neglecting all other terms and imposing 
oj = uji + Atui we obtain 


oo oo oo 


Ezi = 


lUlo 


cz 


dlx / dly / (iUp(U,/, ca) exp[i<ho] exp m 

Ui 


. Awi 


X E 

m=—oo 


— oo —oo —oo 

oo r r TF ^ 


- lJm(u)J2-2m(v) exp 


i[2 — 2m]7r 
2 


-Jm{u)J-2m{v)exp 


i[—2m]7r 


2 m(^) j 

+ 9yJm{u)Ji-2m{v)exp 


i[l — 2m]7r ] 

2 j 

f i[l — 2m]7r 


(34) 


For any value of K and 9^ much smaller than I/ 72 , v is a small parameter. 
This allows one further mathematical step, that is the expansion of the Bessels 
functions Jq{v) ~ in Eq. (27). If one retains only the smallest indexes, one 
obtains the usual expression for the first harmonic of the undulator held under 
the resonance approximation, which is linearly polarized along the eT, direction. 
This can be done by keeping the current term with m = 1 (that is the hrst of 
the two terms in the e7 direction), and neglecting everything else. However, 
here we want to investigate corrections to the held along the Cy direction. 
With this in mind, we also retain gradient terms corresponding to m = 0 and 
m = — 1, which correspond to the largest corrections along the Cy Eq. (34) 
then yields: 


1 ,, = ^^ 


cz 


K ^ 2777 


dL 


dly / dz' exp [i<ho] 


— CXD —CXD 


Xp(U, /, oj) exp[iUU] exp iyf + 9'^ ky^z' 


(35) 
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where we have dehned 


(2(2 + ir2)) '^^{^2{2 + K^^ ' 

^ (2(2+ 7^2)) + "^1 (2(2+ 7^2)) ’ 


and 


$ = 


C-^lo 


-{OJa: + 6*2;/^) + -(^a: + Oy] 


Eq. (35) can be also written as: 


(36) 

(37) 


(38) 


~ Wlo 
h±i = — exp 
cz 




00 CXD 00 


X J J J dz'exp 

—00 

X exp 


K ^ 2K'y 

^AjjCa; + 2 _|_ J^2 ^JjdxdyCy 


— i-— {dxlx + dyly) 


— 00 —00 —00 


p[z ,l,uj) exp[iCz] , 


(39) 


Note that usually computer codes present the product p(x', /, u) exp[iCx'] com¬ 
bined in a single quantity tipically known as the complex amplitude of the 
electron beam modulation with respect to the phase tjj = k^z' + {u/c)z' — ut. 
Regarding such product as a given function allows one not to bother about a 
particular presentation of the beam modulation. Eq. (35) or, equivalently, Eq. 
(39) are our most general result, and are valid independently on the model 
chosen for the current density and the modulation and can be used together 
with EEL simulation codes for detailed calculations of the evolution of the 
vertically polarization contribution to the EEL radiation. 


2.1 Discussion: other contributions to the vertically polarized field 


Inspection of Eq. (35) shows that the vertical polarization component depends 
on angles and roughly scales, aside for the R'-dependent correction factor Rjj 
of order unity, as 2Kfizdxdy Therefore, the ratio between vertical and horizon¬ 
tal polarization components scales as 2fi‘^6x0y. It is useful to underline that the 
vertical-polarization correction was calculated on the basis of the leading reso¬ 
nant contribution to the vertically polarized held. A natural question pertains 
the magnitude of the non-resonant corrections, which -roughly speaking- are 
expected to be suppressed due to fast oscillations in the exponential function 
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in the integrand of Eq. (27). The leading non-resonant correction term is found 
in Eq. (27) setting = ±1 and choosing n = 0, which is justified because 
w 1. The oscillating integrand in dz' suppresses the correction term of a 
factor proportional to the inverse number of undulator periods involved in the 
EEL process, that is those in a held gain length, iV“^, and the correction term 
is therefore proportional to 9y/N^ aside, again, for the i7-dependent correction 
factor Bjj of order unity. Therefore, altogether, the ratio between the lead¬ 
ing non-resonant correction term and the horizontal polarization component 
scales as ^z9y/Nw It follows from the ratio with the horizontally polarized 
main component of the radiation that the non-resonant correction terms can 
be neglected whenever 3> 1. This means, for example, that on-axis 

at 6'a; = 0 the hrst non-resonant term is larger than the hrst resonant term, 
and seems in contradiction with our choice to neglect non-resonant terms. 
However, here we are interested in the different polarization contributions to 
the angle-integrated powers, and not to the intensity at a hxed angle. As al¬ 
ready discussed, the presence of the EEL process limits the detuning to values 
C ky, and the angles of interest up to ~ l/{Ny,^l) l/7z- From 

Eq. (35) follows that the contribution from the leading resonant term is of 
order d(j) d6 6^^ cos^(0) sin^(0) = 7 r/( 67 ^A'^), while from Eq. (27) 

one obtains that the contribution of the leading non-resonant term is of or¬ 
der 72 Jq"" dcf)dO 9^ sin^(0) = Therefore, the non-resonant 

contribution to the angle-integrated polarization correction is fully negligible, 
because the ratio between resonant and non-resonant contributions scales as 
Nyj 3> 1. Similar reasoning can be made, starting from Eq. (19), for the case 
of any other non-resonant contribution arising from any part of the electron 
trajectory, not only the undulator. The specihc angular dependence will vary 
from the kind of contribution but in general one will obtain a held scaling 
as 1/Nyj because it is non-resonant. The transverse distribution extends over 
angles of order ^/W^ times larger than 6'max for the same reason. Therefore, 
one can generalize the argument above, saying that, as long as one consider 
an angular acceptance of order ^max, the power contribution due to any non¬ 
resonant terms will scale as l/Aj); and will be negligible. 

Summing up, in the case of an EEL, due to the presence of a maximum angle 
6 *max related with the self-consistent process, the angle-integrated correction to 
the power from the horizontally polarized radiation component only includes 
the leading resonant term, and Eq. (35) can always be used to calculate such 
correction at the hrst harmonic. 


3 Physical situations treatable analytically 


In the previous section we calculated the two main diherent polarization con¬ 
tributions to the radiation from an EEL with a planar undulator. One orig- 
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inates from the horizontal component of the electric held and the other 
from the vertical component Ey, respectively the a-mode and the 7r-mode. 

We now restrict our attention to the steady-state model of an FEL ampliher. 
As we will see in the hrst part of this section, the simplicity of the steady state 
model offers the opportunity for an almost complete analytical description of 
the polarization characteristics of an FEL ampliher in the high-gain linear 
mode of operation. The second part of the section is devoted, instead, to an 
analytical solution for the case of a constant electron density modulation. 
This model is important from methodological point of view, since it allows 
one to imitate the polarization characteristics of FEL radiation at saturation. 
These analytical solutions can serve as a reliable basis for the development of 
numerical methods. 

Because of the steady state assumption we restrict our attention to one sin¬ 
gle frequency. This means that, in the time domain, the electric held en¬ 
velope E±i must correspond to a real electric held at a certain frequency 
u = oJioil + Ckyj) given by E{z,f±,t) = Tj_i( 2 :, rj_) exp[ia;( 2 :/c — t)] + C.C., 
where the symbol C.C. indicates complex conjugation. In the following we 
will be interested in integrating the angular power distribution of radiation to 
obtain the power fractions into the two modes of polarization. The power for 
the (7- and 7r-polarization components of the hrst harmonic radiation can be 
found remembering that the Poynting vector S = c/{4:7i)E x B represents the 
energy how per unit time and per unit area. The average power carried by an 
electromagnetic signal per unit area can be found by averaging over the signal 
duration. For a monochromatic signal, we further need to take a limit for an 
inhnitely long signal. Distinguishing between the two polarization modes this 
procedure amounts to 


dIF((7,7r) 

dS 


An T^oo T 


T/2 


dt 




(40) 


-T/2 

where we used the fact that in the paraxial approximation and in cgs units 
\E\ = \B\. From our previous discussion we can write 


d^{x,y) 


= 4 


S 




cos^[a;iot + (j){z) 


Substitution in Eq. (40) yields 


( 41 ) 


dIF((7,7r) 

dS 


n 


S 


T/2 

2 1 /* 

±i{x,y){z,f±) Jim — / dtcos‘^[uiot +(piz)] 

T^oo J J 


-TI2 
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c 


2 


( 42 ) 


Finally, integrating over dS^ one obtains the average power for the two frac¬ 
tions 

oo oo 

= ^ J dx J dy\S±i(^,y){z,x,y)\^ , (43) 

— OO —OO 

In order to find S±i(^x,y) we note that the Fonrier transform of onr real monochro¬ 
matic held is given by 


E±{z,f±,u) 

= 27rexp (^i—z ) \£i_i{z, rj_)5(a; — a)) -|- £]_i{z, + <x) 


(44) 


where the last eqnality follows from direct application of the Fourier transform 
shift theorem, since the dependence of the monochromatic held on time is 
zjc — t. By dehnition of held envelope one has thus 


E±i{z,r^,uj) = 271 £±i{z,f±)6{u - ca) -h £li{z,r±)S{u u) 


(45) 


Passing to the complex notation for the held in the time domain (i.e. dropping 
the complex conjugate term), one can simplify the previous equation keeping 
only the positive frequency u and write 


E±i{z,r±,u) = 27i£±i{z,r±)S{u — u) . (46) 

Summing up, in order to calculate IF(o-, 7 r) we should make use of Eq. (43). In 
order to do so, we need an expression for £±i(x,y), which can be found in terms 

of Ej_i with the help of Eq. (46). Finally, one needs to calculate E±i, which 
can be done using Eq. (39). 

Let us therefore turn to the calculation of E±i in Eq. (39), which requires spec- 
ihcation of the bunching p. Formally, the one-dimensional steady state theory 
of FEL amplihers deals with the amplihcation of a plane electromagnetic wave 
by an inhnitely wide and inhnitely long electron beam. The approximation of 
an inhnitely long electron beam is acceptable when one considers a uniformly 
dense beam that is much longer than the slippage in a gain length. However, 
even in this case the electron beam and the electromagnetic wave have h- 
nite transverse dimensions, and dihraction ehects always take place. For the 
one-dimensional model to be applicable dihraction losses must be negligible. 
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In practice, such assumption is valid for XFELs operating in the hard X-ray 
wavelength range. In this case, the asymptote for the held growth rate of the 
fundamental TEMqo mode can be found from the cubic eigenvalue equation of 
the one-dimensional model of EEL ampliher. We also assume that the trans¬ 
verse current density of the electron beam can be modeled as a Gaussian. 

With these assumptions, we may write the charge density in the 

time domain as 


p{z,f^,t) =—p^{r^-roi)f . (47) 

W V VoJ 

The quantity pj_ has the meaning of transverse charge density distribution, 
while / is the temporal electron number density distribution. If the bunch is 
modulated at a given frequency a), one has 


/ 






< a(z )exp 

iu ( 


1 

L V 

VoJi 


Taal) (48) 


where /o is the temporal charge density distribution of an unmodulated bunch, 
assumed to be constant along the undulator, while a is the modulation or 
bunching amplitude. Note that in the general case, the bunching amplitude is 
a complex function of the position in the undulator ; 2 . 


We start the analysis for the case of a stepped prohle electron pulse of hnite 
duration T and then we go over to the limit of an inhnitely long pulse. The 
beam current can be written in the form: 


I 


{-e)N 

T 


Ht 




< a(z) exp 

iu ( 


1 

L V 

Vo/i 


+ C.C. 


(49) 


where HT{t) = 1 in the range (—T/2,T/2) and zero otherwise, N is the 
number of electrons in the pulse and {—e)N/T = —Jq is the average beam 
current. Comparison with Eq. (48) yields 


/ 




I a ( 2 ;) exp 


iu 



+ c.c.|) 


(50) 


Assuming a Gaussian transverse charge density distribution of the electron 
beam with rms size a given by 


pi(f±) 


(-e)N 

2TTa‘^ 


exp 



(51) 


17 



we thus obtain 


p{z,r±,t) 


1 lo 

Vz 27rcr^ 


exp 



Ht 



X 



^a{z) exp 


iuj 



+ c.cX'\ 


(52) 


Let us study the asymptote for an inhnitely long electron beam, that is the 
limit for T —)■ cx), N —)■ oo and {—e)N/T = constant = —Jq. One needs to 
calculate the temporal Fourier transform of p. Fourier-transforming Eq. (52) 
and using Eq. (9) we write the expression for the slowly-varying amplitude 


p{z,f±,u}) = ‘^^^^^-^2na{z)6{uj — uo) , 


where we dehned the current density 


jo{r±) = 


2Tra‘^ 


exp 


2a2 


(53) 


(54) 


and where we dropped the term in (5(a; -|- co) passing to complex notation, as 
done before with the held. 


3.1 High-gain linear regime 


We hrst model the case of an FEL ampliher in the high-gain linear regime. 
We proceed approximating the detuning parameter C as constant along the 
undulator. Let us restrict, for simplicity, to the case of perfect resonance for 
C = 0. This means that from now on a; = Uio- The high-gain asymptote of 
the one-dimensional steady-state theory of FEL amplihers yields 


a{z) = afexp[{V3-\-i)z/{2Lg)] , (55) 

where we set the exit of the undulator (in the linear regime) at z = 0 and 
ttf = constant is the modulation level at z = 0. Here Lg is the held gain 
length. The number of undulator periods in the held gain length Lg is just 
Nyj = (47rpiD)“^, where the FEL parameter pio [9] is related to the problem 
parameters through 


PlD 


47r 




1/3 


1 


(56) 
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where I a = rne(? je ~ 17 kA is the Alfven current. If we now integrate Eq. 
(39) from = —oo to 2 ; = 0 we get 


g 2TiafUio 

Eai = - i -exp 


c^z 


OO 00 0 




K ^ 27^7 


X J dlx J dly J dz'exp 


{ej, + Oyly) 


— OO —00 —00 


X exp 


exp 


(a/3 + i)z' 


2 L. 


^(ca-cjio), (57) 


where we substituted Vz with c in Eq. (53), based on the fact that I/ 7 I 1. 


If we now substitute Eq. (54) in Eq. (57) and we perform all integrations, Eq. 
(46) yields 


-7 2/oOjcaioAg 

£ai = -7^-^ exp 


c^z 




X 


exp 


■'^Mo {01 + »?) /( 2 c" 


K ^ 27^7 

^AjjCa; + 2 _|_ 


(v/i + i)c + iLgUio {O"^ + 6y 


(58) 


The next step is to calculate the two power fractions corresponding to the 
a and tt polarization modes using Eq. (43). It is convenient to present the 
expressions for and HA as a function of the Fresnel number TV, dehned as 


N 


CLin 


(59) 


Since we integrated along z' from — cx) to 0, we did not need to define explicitly 
the length of the undulator. However, in order to compare the for HA and HA 
(and more precisely their ratio HA/HA) with the results obtained above, we 
still need to present the power as a function of the Fresnel number N, dehned 
exactly as before. 

Explicit calculations yield 

VH^.y VBjjViDGAiV) J 

where 
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G, 


-i.N) = ^ exp[(l - iV3)N] {tt + Trexp [2iV3N] 

2iy/3N] Ei (iV(-l - iVs)) + iEi {iN{i + ^3)) } 


-i exp 


( 61 ) 


— exp[(l + iy/3)N] \n — iEi ((—1 — iv^)iV^j 

— (Si + exp[(l — iy/3)N] \n + iEi [i{i + | (62) 


and 


with Wh = meC^'jIo /e is the total power of the electron beam. 

The ratio between the fractions radiated in the two modes of polarization is 
therefore conveniently expressed as a function of three separate factors: 


= f(K)g(N^)w{N) 

with 


(64) 


fiK)- 

9{Pid) 
w{N) : 


till 

pId 

'G^{N) 

G.{N) 


(65) 


The hrst factor, f{K), is only a function of the undulator K parameter and 
is plotted in Fig. 1. As it can be seen from Fig. 1 it always remains between 
about unity and 3, for any value of K. The second factor, g^pin) scales as the 
inverse number of undulator periods squared, and is a signature of the fact 
that the gradient term in the equation for the electric held scales as the inverse 
number of undulator periods. The third factor, w{N), is only a function of the 
diffraction parameter that is, once the wavelength and the undulator length 
are hxed, a function of the electron beam size only. It is also plotted in Fig. 1. 
It is unity for values of the diffraction parameter around unity, but it quickly 
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Fig. 1. Illustration of the behavior of f{K) (left) and 'w{N) (right). 

decreases for larger values of N. Clearly the power fraction radiated in the 
71 mode increases drastically with the photon energy, partly due to a larger 
number of undulator period per held gain-length, but mainly because of a 
larger diffraction parameter. 

For the sake of exemplihcation we apply the model in the linear regime dis¬ 
cussed in this paragraph to one particular case, calculated from start-to-end 
simulations for the SASEl and SASE2 line of the European XFEL. We con¬ 
sider a 250 pC electron beam at a photon energy of about 9 keV. The electron 
energy is 17.5 GeV and K ~ 3.6. The peak current is about 5 kA, and the rms 
sizes of the electron beam in the horizontal and vertical directions are about 
ax — 15 yum and ay ~ 18 /rm respectively. For our purposes of exemplihca- 
tioup] we consider a round beam with a = 16/im. The peak current density 
can then be estimated as Jo/(27r(T^). Finally, the undulator period is = 40 
mm. From these numbers we obtain the parameter piD ~ 8 ■ 10“^. Plugging 
these numbers in Eq. (72) and remembering the dehnition in Eq. (59)we ob¬ 
tain f{K) ~ 2.5, g{piD) — 6.4 ■ 10“’^, AI ~ 3 and w{N) ~ 0.072, so that the 
overall ratio Wt^/W^ ~ 1.13 ■ 10“^. 


3.2 Constant density modulation 

In analogy with the previous paragraph, we now proceed to study the case of a 
constant density modulation along an undulator of hxed length imitating 
the behavior of an FEE at saturation. We can still set C{z) = 0. At variance 
with the previous model we now write 



( 66 ) 


^ If the electron beam is not round, one can easily modify the model for the electron 
beam distribution in this paragraph. 
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Here a/ = const a constant modulation level, = 1 for 2 ; in the range 

{—Ly^/2, L^/2) and zero otherwise, with the undulator length, and jo is 
dehned as in Eq. (54). In this case, Eq. (39) can be written as 


= 2nafUio 
E±i = -^-exp 


c^z 


.^1 




^ A ^ 
T^AjjCj; + 


2K'j 


X 


CXJ CXJ CXJ 

J dlx J dly J dz' exp 


—i- 


27 

.CJlo 


2 + K^ 

i^dxlx “h 


^JjdxOyCy 


—00 —00 

'.U)l 


X exp 


^ do (f) HlAz')5{uj - uu 


(67) 


The integral in dl amounts to the spatial Fourier transform of jo (jj ■ We 
proceed similarly as in the previous paragraph and obtain 


Sn = 


loCl ^UJ\oL^i 
c^z 


exp 




X 


K ^ 2K-i 




xsinc 


L-u)(jJ\o 

4c 


[el + 0. 


exp 


2c2 


(el + 0 


( 68 ) 


The hnal step, as in the previous paragraph, consists is the calculation of 
the angle-integrated hrst harmonic power. As before the power for the x- and 
y-polarization components of the hrst harmonic radiation are given by Eq. 
(43). 


Both integrals for the horizontal polarization component and the vertically 
polarized correction can be calculated analytically by exploiting the cylindrical 
symmetry of the model. One is then left with 


[wj °I^Bjj2(47riV^)-iF^(iV)J ’ 


where 


Fo{N) = arctan 



+ iVln 


/ 4iV2 \ 
V4iV2 + 1 J ’ 


(70) 


F^N) 


1 

iV(l + 4iV2)’ 


( 71 ) 
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Fig. 2. Illustration of the behavior of h{N). 
and where parameters N and Wo are given by Eq. (59) and Eq. (63). 

Similarly as before, the ratio between the two fractions radiated into the two 
modes of polarization is conveniently expressed as a fnnction of three separate 
factors: 

= f(K)g(N^)k(N) (72) 

with 


fiK) = 
9{N^) = 
h{N) = 


1 

'F^{N) 
FAN) ■ 


(73) 


The fnnction / has been defined in the previons paragraph. Concerning the 
second factor g, we have an expression which is similar to that in Eq. (65). 
The only difference is that here we replaced pm with (47riV.u,)“^, with N^, 
the nnmber of nndnlator periods in the nndnlator. The nnmber of nndnlator 
periods in a held gain length is jnst = (47rpiD)“^, and therefore the second 
factor in Eq. (72) jnst amonnts to for an nndnlator length = Lg, Lg 

being, as before, the held gain length. By setting the nndnlator length eqnal 
to the held gain length the two models can be directly compared by stndying 
w{N) as dehned in Eq. (65) and h{N) dehned in Eq. (73). We plot h{N) 
explicitly in Fig. 3.2. As one can see it dihers from Fig. 1, dne to the diherent 
model nsed. 


Considering the same example made in the previons paragraph we hnd again 
f{K) ~ 2.5, gipio) — 6.4 ■ 10“’^, ~ 3. Plngging the valne for N into Eq. 

(73) we obtain h{N) ~ 0.097, so that the overall ratio W-^/W^ ^ 1.5 ■ lO”"^. 
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3.3 Polarization characteristics for the second harmonic 


So far we have discussed the values of the two polarization components a and 
TT. We have shown that, typically, in the case of an XFEL with a horizontal 
planar undulator, only less that one part in a million of the total power at 
the hrst harmonic is polarized in the vertical direction. For some experiments 
even such small fraction of the vr mode is of importance. 

In addition to the radiation at the hrst harmonic, however, there are several 
other background contributions to the vertically polarized power fraction from 
higher harmonic radiation. It is useful to be able to estimate what is the typical 
fraction of the vr mode component at higher harmonics. There is an important 
distinction to be made between odd harmonics and even harmonics. The vr- 
mode contribution from all odd harmonics {uj = Scuio, oj = dwio, etc.) can be 
completely disregarded, while in practice even harmonics {uj = uj = dwio, 
etc.) can still be an important source of vertically polarized radiation. Since 
the total power contained in each harmonic is proportional to the square of 
the bunching amplitude, only the contribution due to the second harmonic is 
of practical interest. 

The contribution of the second harmonic can be calculated using results in 
[To] . We proceed to study the case of a constant density modulation. The 
ratio between the fraction of total power radiated into tt mode at a; = 2ujio 
and radiated into the a-mode at a; = ujio can be expressed as a function of 
separate factors [TO] : 

^ 1 2 + K^ al F2^(N) 

47rN^ a)A^jjF,{N) ’ ^ ^ 

where N and are given, respective!}]^ by Eq. (59) and Eq. (70), Ajj is 
given by Eq .(36), 02 and a/ are the amplitudes of the beam modulation at 
u = 2ujio and uj = ujio respectively. 


and 

(1 + I^) (76) 

^ Note that N in [10] is defined twice larger with respect to what is reported here. 
In fact, in [TO] all results refer to the Fresnel number for the second harmonic. 
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Fig. 3. Illustration of the behavior of F 2 {N)/Fa-{N). 

We plot F 2 {N)/F^{N) explicitly in Fig. 3. Considering the same example made 
above we find 


VFi. 


10 -"^ 

aj 


(77) 


The contribution from the even harmonics scales therefore as can 

be completely disregarded when the XFEL operates in linear regime. At satu¬ 
ration, the contribution from the second harmonic can be comparable with the 
first harmonic in the case when X-ray optics harmonic separation is absent. 


4 Conclusions 


One attractive feature of radiation from X-ray Free-Electron Lasers (XFELs) is 
its high degree of polarization. This paper shows that for an FEL with a planar 
undulator with the electron motion on the horizontal plane, the horizontally 
polarized component of radiation from greatly dominates the photon beam 
characteristics and only less than one part in a million of the total intensity 
is polarized in the vertical plane. This feature that can be useful in different 
experimental situations. When describing physical principles, analytical de¬ 
scriptions are always important: they allow for a proper understanding of the 
principles, in our case those of FEL physics, and also for testing numerical 
simulation codes. From this point of view, a SASE XFEL is a rather com¬ 
plicated subject. In fact, in all generality, its radiation can be represented as 
a non-st at ionary random process whose analytical description is complicated 
by the fact that the electron bunch combines both the features of the input 
signal and of the active medium with time-dependent parameters. Approxi¬ 
mations are therefore needed. In particular, it is important to find a model 
allowing for an analytical description without loss of essential information 
about the features of the FEL process. A model satisfying these conditions 
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is that of a stepped profile electron bunch, together with the application of 
self-seeding scheme for narrowing the radiation bandwidth. In the framework 
of this model it becomes possible to describe the polarization properties of the 
radiation from an XFEL in a fully analytical way. Using Maxwell’s equations 
one can write an explicit expression for calculating the electric held with given 
charge and current sources. Up to now, in all FEL codes the contribution of 
the charge term is assumed to be negligibly small. However, in our case of 
interest, the charge term is the only responsible for the vertically polarized 
radiation component. Our analytical results, in particular those for the high- 
gain linear regime, are therefore expected to serve as a primary standard for 
testing future FEL codes including the charge term as electromagnetic source. 
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